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^ Abstract. A set of gauge invariants are identified for the gauge theory of quantum 

Cd anholonomies, which comprise both the Berry phase and an exotic anholonomy in 

eigenspaces. We examine these invariants for hierarchical families of quantum circuits 
whose qubit size can be arbitrarily large. It is also found that a hierarchical family of 
quantum circuits generally involves an NP-complete problem. 
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1. Introduction 

An adiabatic time evolution of a quantum system along a closed loop may induce a 
nontrivial change. The best known among them is the phase anholonomy, which is also 
referred to as the Berry phase: a stationary state acquires an extra phase of geometrical 
origin as a result of an adiabatic cycle [1]. It has been recognized that the eigenvalue 
and eigenspace of a bound state can be also changed by an adiabatic cycle [2J: When we 
keep track of an eigenvalue along a cycle, discrepancies may appear between the initial 
and finial eigenvalues. Such a change, which is called an eigenvalue anholonomy, implies 
the interchange of the eigenvalues, which is in fact allowed because of the periodicity of 
the spectrum for the cycle. Furthermore, due to the correspondence between eigenvalues 
and eigenspaces, the eigenvalue anholonomy implies the eigenspace anholonomy, where 
an adiabatic cycle delivers one stationary state into another stationary state. 

The earliest example of the eigenvalue and eigenspace anholonomies is found in a 
family of Hamiltonians with generalized point interactions [3]. Later, it is recognized 
that there is another examples of anholonomy in Hamiltonian systems that involve 
crossings of eigenenergies [H [3]. Further examples are found recently in periodically 
driven systems and quantum circuits, where a stationary state is an eigenvector of a 
unitary operator that describes the evolution for a unit time interval |6l |71 E]- The 
spectral degeneracies are known to make both the phase anholonomy |9] and the 
eigenspace anholonomies more interesting. However, throughout this paper, we will 
focus on the case where no spectral degeneracy exists along adiabatic cycles. 

We have recently developed a unified gauge theoretic formalism of eigenspace 
anholonomy [TU], extending the Fujikawa's approach to the phase anholonomy [TT]. 
The key quantity for the formalism is the holonomy matrix M{C), whose elements are 
given by overlapping integrals between stationary states and their adiabatic transports 
along a cyclic path C. Because M{C) is gauge covariant [lO], only a part of M{C) is 
observable [T2] . 

In this paper, we show how we extract gauge invariants from M{C). Let us suppose 
that all stationary states in A^-dimensional Hilbert space are connected by the A^- 
repetition of the adiabatic cycle. It can be shown that M{C), under an appropriate 
choice of the gauge, is given in terms of two gauge invariants &{C) and 7(C): 

M(C) = 6(C)e'^(^)/^, (1) 

where &{C) is a permutation matrix that describes the interchange of eigenspaces 
induced by the adiabatic cycle C, and 7(C) is the Berry phase gained through the 
iV-repetition of the cycle C. In general, M{C) is decomposed into the blocks, each of 
which consists of a permutation matrix and a geometric phase factor as shown in ([I]). 

Also, in this paper, we provide examples of systems that exhibit various types 
of &{C) by extending a recursive construction of hierarchical quantum circuits [8j. 
Due to its topological nature, 6(C) remains unchanged under small perturbations, and 
the eigenspace and eigenvalue anholonomy with various &{C) are robust against small 
imperfections inevitable in experimental implementation. 
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The plan of this paper is the following. We first lay out the concept of the eigenspace 
and eigenvalue anholonomies with the minimal model, a quantum circuit on a qubit, in 
Sectional In Sectional we revisit our gauge theoretic formalism [10] to reveal the hidden 
gauge invariants &{C) and 7(C) residing in M{C). In the following sections, we attempt 
to extend the recursive construction of quantum circuits [8] in order to obtain expanded 
instances of families of systems with variety of ©(C). In Section 111 we introduce novel 
families of quantum circuits. The eigenvalue problem of these circuits are also solved. In 
Section [5} we examine the eigenvalue anholonomy of our hierarchical models, and obtain 
(5(C) for each family of quantum circuits. In Section [61 we examine the eigenspace 
anholonomy of the hierarchical models to obtain 7(C) hidden in M{C). In Section (tI 
several examples are shown. We examine the relationship between our model and the 
adiabatic quantum computation [13] along the eigenangle [H] in Section [s] In Section [9} 
we discuss the relationship between our result and a topological characterization of the 
winding of quasienergies for periodically driven systems [15] . A summary of this paper 



is found in Section 10 Several technical details are provided in an Appendix. 



2. Anholonomies in quantum circuits 

We lay out the eigenvalue and eigenspace anholonomies and associated gauge invariants 
using a family of quantum circuits on a qubit. These quantum circuits provide the 
simplest case for our analysis shown in latter sections. 

We introduce a quantum circuit on a qubit, or, equivalently a quantum map for 
spin-^ [ini[IS]: 

uiX,p)^exp{i{p~l)\{l-\y){y\)+i\\y){y\}Z, (2) 

where p is an integer. Using orthonormal vectors |0) and |1) of the qubit, we define 
\y) = (|0) — i\l))/\/2 and Z = |0)(0| — |1)(1|. Because ^(A,^) is 27r-periodic in A, we 
will examine the eigenvalue and eigenspace anholonomies for the cycle C, for which A 
is increased from to 2tt. We denote an eigenvalue of the unitary operator u{X,p) as 
exp{i^(n; A,p)}, where an eigenangle 9{n] X,p) satisfies 

9{n;\,p) =n7r+ -X, (3) 

for n = 0, 1. The corresponding eigenvectors are the following: 

|0(A,p)) = |0)cos^^^^ + |l>sm^^^, (4a) 

|l(A,rt) = |l)coB^-|0)=in<i^. (46) 

First, we examine an anholonomy in the eigenangle 6{n]X,p). As A is increased 
from to 2tt along the closed cycle C, we obtain the following change 

e{n; X + 27T,p) = e{s{n; C); X,p) + 27rr(n; C), (5) 

where we introduce integers 

s{n;C) = {n + p) mod 2, r{n;C) = [{n + p)/2] , (6) 
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and [x] is the maximuni integer less than x. Hence, when p is even, n-th eigenangle is 
transported to itself by the cycle C. On the other hand, odd p implies the presence of 
an anholononiy in the eigenangles, i.e., 6{n; A + 2iT,p) = 6{n; X,p) + 27m, where = 1 
and 1 = 0. 

The eigenangle anholonomy implies an anholononiy in eigenspaces. Namely, 
the adiabatic transport of an eigenspace along the cycle C induces the change in 
eigenvectors. Suppose that each eigenvector \n{X,p)) satisfies the parallel transport 
condition within the eigenspace [17], i.e., 

{n{X,p)\[dx\n{X,p))] = Q. (7) 



This condition is, in fact, satisfied by the eigenvectors (4a) and (46). The adiabatic cycle 
C accordingly transports \n{0,p)) to \n{27i,p)). Because of the correspondence between 
the eigenangle and the eigenspace, |n(27r,p)) must agree with \s{n;C){0,p)) up to a 
phase factor. When p is odd, |n(27r,p)) is orthogonal to \n{0,p)) because of s{n) ^ n. 
Note that the change of the eigenvector occurs in spite of the absence of the spectral 
degeneracy along the path C. Hence this anholonomy is different from Wilczek-Zee's 
phase anholonomy. Instead, an anholonomy appears in the eigenspace. 

The change in the eigenvectors is characterized by the overlapping integrals between 
the initial and the finial eigenvector 

{M(C)}„,,„ = (n'(0,p)|n(27r,p)), (8) 

which we call a holonomy matrix. This involves the interchange of eigenspaces and the 
accumulation of phases. Although we have assumed the parallel transport condition for 
each eigenspace in (|8]), M{C) is gauge covariant with respect to the gauge transformation 
for |n(A,p)), as will be explained in Section [s] [TO]. Due to the gauge covariance, only a 
part of M{C) is observable [12] . 

We now look for the gauge invariants, i.e. observables, of the quantum anholonomies 
for the adiabatic cycle C . They can be identified easily by investigating M{C) closely. 
One is the permutation matrix ©(C): 

{6(C)}„/,n = 5n',s(n;C), (9) 

which reflects the interchange of eigenvalues and eigenspaces. Other gauge invariants 
are the geometric phases. When M{C) is diagonal, i.e., the eigenspace and eigenvalue 
anholonomies are absent, the n-th diagonal element of M[C) is the Berry phase for 
the n-th eigenspace. On the other hand, if the eigenspace and eigenvalue anholonomies 
occur, M{C) contains the off-diagonal geometric phases, which have been introduced 
to quantify the phase anholonomy for a noncyclic-path [18] under the presence of the 
eigenspace interchange [T^ . 

Let us look at the holonomy matrix and the gauge invariants in our model ([2]). 
Since eigenvectors (4a) and [Ah] satisfy the parallel transport condition (JTJ), we obtain 



M(C) = lcos^^^-.Fsin^^^, (10) 



Gauge invariants of eigenspace and eigenvalue anholonomies 5 

where 1 is the identity matrix (which will be omitted below) and 

-(!:)■ '"' 

Because the integer p is a crucial parameter for the anholonomies as shown above, 
it is better to classify the cases depending on whether p is even and odd to investigate 
the gauge invariants. When p is even, we obtain s{n; C) = n for n = or 1. Namely 
the anholonomies apparently disappear, i.e.. 




e(c) = . : • (12) 



The permutation described by ©(C) consists of two disconnected cycles each of which 
has a unit length. The holonomy matrix takes block-diagonal form, whose diagonal 
blocks are two 1x1 matrices, i.e., 

where a{n]C) = (— 1)(2-p)/2 jg ^^j^g Berry phase factor acquired through the adiabatic 
cycle C for the n-th eigenspace. 

On the other hand, odd p implies s(0; C) = 1 and s(l; C) = 0. Hence the eigenvalue 
and eigenspace anholonomies appear explicitly, i.e.. 



e(C^) = K n h (14) 




which is the permutation matrix of the cycle whose length is 2. The holonomy matrix 
can be regarded as a product of ©(C) and a diagonal unitary matrix, i.e., 

MiO ^ S(C) ( '<^;f^ ^,^«^j ) . (15) 

where we set a{n; C) = (— 1)"+(i^p)/2. Each a{n; C) has no geometric significance since 
it depends on the choice of the gauge. We can thus construct a gauge invariant quantity 
o"(0; C)cr(l; C) = — 1. This is the Manini-Pistolesi's off-diagonal geometric phase factor, 
which is an extension of Berry's phase factor for the case that the permutation of 
eigenvectors takes place along an adiabatic evolution ^U\. Once we slightly change the 
phase of eigenvectors as \0{X,p)) = i\0{X,p)) and |l(A,p)) = |l(A,p)), the holonomy 
matrix in the corresponding gauge takes a canonical form, i.e., 

M(C) = ©(C)e*^(^)/^ (16) 

where 7(C) = np agrees with the Manini-Pistolesi off-diagonal geometric phase, up to 
a modulo 27r. 

So far, we have examined specific families of quantum circuits, namely u{X,p) (pi) 
with \y) = (|0) — i|l))/A/2 and Z = |0)(0| — |1)(1|. This analysis can be applied to the 
whole family of u{X,p) with arbitrary \y) and Z. For the path Cp of the whole family 
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u{X,p), Cp still has topological invariant, namely either (12) or (14). Although we may 
deform the closed path Cp by varying \y) and Z in u{X, p), the topological invariant &{C) 
remains intact unless the deformed path encounters a spectral degeneracy. Hence, for 
every path C obtained by deforming Cp with even p, &{C) is the unit matrix. With 
odd p, 6(C) is always the permutation matrix of the cycle whose length is 2. Hence 
the closed paths are divided into two families related to the value of ©(C), i.e. the unit 
or the permutation matrix except the paths that come across spectral degeneracies. 

3. Gauge invariants of the quantum anholonomy 

In this section, we first review our unified gauge theoretical treatment of quantum 
anholonomy that is applied both to the Berry phase and the eigenspace anholonomy [TOl 
[TB] . This offers a gauge covariant expression of the holonomy matrix M{C) using a non- 
Abelian gauge connection. Second, we reveal gauge invariants in M{C) using several 
choices of the gauge. In particular, we will show that Manini-Pistolesi's off-diagonal 
geometric phase in M{C) is the Berry phase of a hypothetical closed path. 

Because the eigenspace anholonomy involves the interchange of multiple 
eigenspaces, we need to treat basis vectors of the relevant eigenspaces simultaneously. 
Our starting point is the sequence of ordered basis vectors 

/(A) = [|0(A)),...,|Ar-l(A))], (17) 

where N is the dimension of the Hilbert space of our system and the basis vectors are 
assumed to be normalized. A non-Abelian gauge connection is then given as 

A{X)^[f{X)]^idJiX), (18) 

where the (n',n)-th element of A{X) is {n' {X)\idx\'n{X)) . The anti-path-ordered integral 
of A{X) along C, 



W{C) = exp ( -i / A{X)dX J , (19) 

describes the parametric change of /(A) along C. Note that W{C) = {/(0)}"''/(27r), 
where and 27r correspond to the initial and final points, respectively, in C. We also 
introduce a diagonal gauge connection 

{AD(A)}„,,„ = 5„,,„{A(A)}„,,„, (20) 

whose 77,-th diagonal element is Mead-Truhlar-Berry's gauge connection [201 [1] for n-th 
eigenspace. If the eigenspace anholonomy is absent, the path-ordered integral of A^{X) 
along C, 

B{C) = exp A f A°(A)dA j , (21) 

contains the Berry phases. 

In general, neither W{C) nor -B(C) is gauge invariant. Indeed, the gauge 
transformation |n(A)) H- e^'^"(^^|n(A)), when the spectral degeneracy is absent, induces 
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the change W{C) ^ G{0yW{C)G{27i) and B{G) ^ G{2nyB{C)G{0), where G{X) is 
the diagonal unitary matrix whose n-th. diagonal element is e"'^"'-'^-' [10] . This means that 

M{G) = W{G)B{G) (22) 

is gauge covariant, i.e., M{G) i— )■ G{OyM{G)G{0). Another proof of the covariance 
of M{G) is obtained through the investigation of adiabatic time evolution of the basis 
vectors directly [TUl [IS] . 

The holonomy matrix M(C) is also expressed as a product of a permutation matrix 
&{G) and a diagonal unitary matrix cr(C), i.e., 

M{G) = e{G)a{G), (23) 

because W{G) involves the interchange of the eigenspaces and B{G) is a diagonal unitary 
matrix. A permutation matrix can be decomposed into cycles, each of which is a cyclic 
permutation [21]. When &{G) is decomposed into J cycles &j{G) (j = 1, . . . , J), M{G) 
takes a block diagonal structure where j-th diagonal component is &j{G)aj{G) with a 
diagonal unitary matrix (Jj{G) under a suitable choice of a sequence of the basis vectors. 
In the following, we will focus on a diagonal block, or equivalently, we assume J = 1, 



i.e., &{G) (23) describes a single cycle. Accordingly, Manini-Pistolesi's off-diagonal 
geometric phase factor is 

e'^(^) = det(T(C). (24) 

Both ©(C) and 7(C) are gauge invariants of the adiabatic cycle G. 

With an appropriate gauge of the basis vectors, M{G) takes a simpler form 

M{G) = 6(C)e'^(^)/^, (25) 

i.e., the off-diagonal geometric phase 7(C) is equally assigned to each eigenspace. This 
tells us that the two gauge invariants &{G) and 7(C) are at the heart of M{G). We 



thus show (25) by using a diagonal unitary matrix Ud that satisfies 

f/]{6(C)a(C)}f/d = 6(C)e'^(^)/^, (26) 



where the left hand side describes a gauge transformation of (23) induced by Ud- First, 
we assume that the basis vectors are arranged so as to satisfy {&{G)}n',n = Sn',n+i mod n- 
This is always possible for a cyclic &{G). Second, we set 

{UdU = 1, 

{f/d}ii = {f/d}oo{^(C)}ooe-'^(^)/^, 

{Ud}22 = {f/d}ii{a(C)}ne-'^(^)/^, (27) 

{Ud}N-l,N-l = {f/d}^-2,^-2{^(C)}iV-2,7V-2e-'^(^)/^. 



Hence it is straightforward to confirm ( |26[ ). 

Now we construct a gauge that precisely assigns &{G) and <j{G) to W{G) and 
B{G), respectively. Let us suppose that such a gauge is chosen. Accordingly W{G) 
describes only the interchange of eigenspaces. On the other hand, the Manini-Pistolesi 
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off-diagonal phase 7(C) is encoded into B{C) because det i?(C) coincides with e 



hiC) 



This gauge offers a way to associate 7(C) with the gauge connection from (21): 



7(C) = E / {M^)}nnd\. (28) 



Note that the argument above does not guarantee that (28) is gauge invariant because 
det B{C) is not gauge invariant in general. Nevertheless we will show that 7(C) is indeed 
gauge invariant below. 

For this we examine the origin of the deviation of W{C) from &{C) in an arbitrary 



gauge. First, let us recall {W{C)}n',n = ('^'(0)|^(27r)) (19). Second, because C is a 
loop in the parameter space, {&{C)}n',n = |{W^(C')}n',nP is either or 1. We introduce 
s{n; C) so that \n{27c)) is parallel to |s(n)(0)), as is done in the previous section. Hence 
(|9| is also valid here. Now we obtain {W{C)}n',n = {6(C)}„/,„,(s(n)(0)|n(27r)). Thus 
we need to impose (s(n)(0)|n(27r)) = 1 using a suitable choice of the phases of the basis 
vectors, to have W{C) = 6(C). 

We introduce a normalized vector |^(A)) to define basis vectors |^n(A)) (n = 
0, . . . , A^ — 1), which satisfy W{C) = &{C). First, for Once we impose smoothness 
on |C(A)), |^(A)) outside < A < 27r is determined. In particular, the projection 
|^(A))((^(A)| is 27rA^-periodic in A. Next, from |^(A)) one can construct |^n(A)) in the 
following way; \U^)) = |e(A)), \Uo){X)) = le(A + 27r)), |e.(.(o))(A)) = |e(A + 47r)), and 
so on. At the final step of this construction, we obtain |^nF(A)) = |^(A + 27r(A^ — 1))), 
where np is defined as s^np) = 0. This imposes (^s(n)(0)|^n(2vr)) = 1 except for n = np. 
The case n = np is related to the single- valuedness of |^(A)) including the phase factor 
in C^, i.e., A^-repetitions along the closed cycle C. Here we simply assume that there 
is a gauge that satisfies |^(A)) is single- valued in C^ . 



Now we return to (28) to obtain an intuitive interpretation of 7(C). In the gauge 



specified by the basis vectors |^„(A)), we have {A(A)}„„ = i(^„(A)|9A|^n(A)). Because 



(28) accumulates the contributions from all n = 0, . . . , A^ — 1 along C, this can be cast 
into an integration along the extended cycle C^: 



7(C) = / (e(A)|[ia,|e(A))]dA. (29) 

Now it is straightforward to see 7(C) is the Berry phase of the adiabatic basis vector 
|^(A)) along C^. Hence it is independent of the gauge of |^(A)), despite our specific 
choice. 

4. A recursive construction of A^-qubit examples 

In the following sections, we extend the quantum circuit u{\,p) ^ to multiple-qubit 
systems. Our aim is to show a various realization of &{C) in many situations. In this 
section, we explain how to extend it and solve the corresponding eigenvalue problem. 

A crucial ingredient of our extension of u{X,p) is the so-called "super-operator" 
Dp[-] labeled by an integer p: 

D,[u]^cy[u]{z®i), (30) 
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y 

T 



[/(A'-i)(;t) -{C/(^-i)(A)}-i+/"v 



Figure 1. The schematic diagram describing the recursion relation (33) of U^^\X). 
The uppermost Hne indicates the iV-th qubit, which plays a role of the control qubit 
of the controlled unitary operations. When the control qubit is one (zero) the unitary 
operation connected to the filled (open) circle is applied to the rest of the qubits. y 
above the circles indicates that the axis of the control-bit is given in the y-direction 
(see (31 1, the definition of Cj^[U]). In the following, y in circuit diagrams is omitted. 



where 

Cy[U] ^ (i - \y){y\) ® U^-' + \y){y\ ® U (31) 

is a generalized controlled-f/ gate. Here the "axis" of the control-bit is given in the 
"^/-direction" . It is worth noting that one can retrieve ii{X,p) from Dp[-] through 

4[f/g(A)]=M(A,p)®iA, (32) 

where Ug{X) = c^^Ia and 1a are the global phase gate and the the identity operator of 
the ancilla, respectively. 

We examine a family of quantum circuits U'^'^\X) for a given sequence of integers 
{Pj}'^=v We define ^/^(A) = u{X,Pi) for A^ = i, and 

f/W(A)^4,[?7(^-i)(A)] (33) 

for A^ > 1 (see Figure II]). If there is no ambiguity, we will omit p^'s. 

We now solve the eigenvalue problem of U^^\X) recursively. Let \nN, ■ ■ ■ , rii, A) and 
9^^\n]\f, ■ ■ ■ ,ni; A) denote an eigenvector and the corresponding eigenangle of U^^\X), 
i.e., 

f/(^)(A)|n;v, ■ ■ ■ , ni; A) = e^''*"'('^-'-'"^^^Viv, ■ ■ ■ , m; A), (34) 

where Uj E {0, 1}. We have already solved the case A^ = 1 in Section [2] (see (|3|, (4a) 
and (46)): I'^i; A) = \ni(X,pi)) and 9^^\ni; X) = 9{ni; X,pi). For A^ > 1, we obtain 
recursion relations for eigenangles and eigenvectors. Suppose we have an eigenvector 
|v|/{^-i)(;^)) and the corresponding eigenangle 0''^~^''(A) of U^^~^\X). Recall that we 
obtain U^^\X) by extending U^^~^\X) so as to incorporate one more qubit into it. Let 
10) denote a state vector for the qubit newly added to the quantum circuit described 
by f/(^"^)(A). A natural candidate of an eigenvector of U^^\X) is |0) (g) |^(^"^)(A)) 
because we have 

f/(^)(A) (10) ® |vl/(^-i)(A))) = HX';p^)m ® |vl/(^-^)(A)), (35) 

where A' = e(^~^)(A). Accordingly, \nN{X',pN)) ® |^(^"^)(A)) and 9{nN;X',pN) with 
un € {0, 1} are an eigenvector and the corresponding eigenangle of U^^\X), respectively. 
This implies the recursion relations 

9^^\nN, . . . ,ni; A) = 9{n^; 9'^''-'\nN.,, ...,n,; X),Pm) (36) 
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and 

\nN, . . . , ni; A) = \nN{e^^~^\nN-i, • • • , ^i; A))) (g) |n7v_i, 
for N > 1. 
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,ni;A), (37) 



Assume that the solutfon of (36) takes the foUowing form: 



e'^^\nN, . . . ,ni; A) = ^ ( mNiuN, . . . ,ni) + (In: 



(3J 



^ I m^(n^, . . . ,ni) + div^ 

A normahzed sfope d^ of the eigenangle and a "principal quantum number" 
mN{nN, ■ ■ ■ ,ni) are defined befow. For the sake of brevity, we will sometimes omit 
the arguments (n^y, • • • ,ni) oim^, so that O^^^n^, . . . ,ni; A) is written as 9^^\m]^] A). 
For A^ = 1, we have ?Tii(ni) = ni and di = pi. For A^ > 1, ( [36) ) provides the recursion 
relations: 



rriN 
Hence we have 

rriN 



2^ ^un +PNmN-i 



and 



dN — PndN- 



d 



N 



N-1 / N 

AT 



fc=i 



The solution of the recursion relation (37) of the eigenvectors is 

N 

\nN,...,ni;\) = (^\nj{9^^'^\mj_i; \),pj)). 



(39) 

(40) 
(41) 

(42) 



We remark on the spectral degeneracy of U^^\\). From (38) and (40), the spectral 
degeneracy is absent if and only if all p2, . . . ,Pn are odd. Note that the degeneracy 
condition is independent of pi and A. Because we focus on the case that there is no 
spectral degeneracy along the path C, we will assume that P2, ■ ■ ■ ,Pn are odd in the 
following. 



5. Eigenangle anholonomy 

We here examine the eigenangle anholonomy of the family of quantum circuits described 



by U^^\X) (33). To achieve this, we keep track of the parametric dependence of the 
eigenangle O^^^njy, . . . , rii; A) along the cycle C, where A is increased by a period 27r. 
The parametric change of the eigenangle can be expressed as 

e^^'^nN, . . . , rii; A + 27r) = e^^\s^^\nN, . . . ,ni; C); A) + 27rr(^)(nAr, . . . ,ni; C), (43) 

where s^^\nN, . . . , Ui, C) is a collection of A^ quantum numbers and r^^^n^, . . . ,ni;C) 
is an integer. In this section, we explain how s^^\nN, ■ ■ ■ ,ni;C) and r^^\niy, . . . ,ni;C) 
are recursively determined. They completely characterize not only the eigenangle 
anholonomy but also the eigenspace anholonomy when the eigenangles are not 
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degenerate. The reason is that eigenangles have one-to-one correspondence to the 
eigenspaces if the spectral degeneracy is absent. On the other hand, because the exphcit 
expressions of these quantities are comphcated in general, we postpone obtaining explicit 
solutions; several examples will be shown in Section [7} 

The recursive structure in U^^\X) determines s^'^^n^, . . . ,ni;C) and r^^^n^, . . . ,ni;C). 
For brevity, we will omit the argument C in the following. Let s, {njy, . . . ,ni) denote 
the quantum number of j-th qubit in s^^^'^niy, . . . ,ni) (1 < j < A^). We note that 
Sj {un, ■ ■ ■ ,ni) is either or 1. For j < N, we obtain 

sf\nN,---,ni) = sf^'^\nN-i,...,ni), (44) 

using the recursive structure in f/*^^)(A). Hence, s^^\nN, ■ ■ ■ ,ni) is determined by 
sf(n„...,ni) (l<j<iV). 

We need to find s)^ {n^, . . . ,ni) and A^\nj\[, . . . ,ni). For A^ = 1, we obtain, 
from the analysis of the single qubit case in Section |2| s[ (^i) = s(ni,pi) and 
A^\ni) = r(ni,pi), where s{n,p) and r{n,p) are defined in (ml). In order to clarify 
the case of A^ > 1, we first examine the change of the principal quantum number rriisf 



by an increment of A from to 27r using (43) and (38) 



mN{,s^^\nN,---,ni)) = mN{nN, ■ ■ ■ ,ni) + dN - 2^r'-^\nN, ■ ■ ■ ,ni). (45) 

Thus we have 

s'-p{nM,...,ni) = s(njv,pjvr(^"^)(rajv_i,...,ni)), (46a) 

r^^\nN, . . . , ni) = r{nN,PNr^^~^KnN-i, • • • , ^i)). (466) 

We note that it is generally difficult to write down the explicit expressions of the solution 
of (46a) and (466). This problem will be dealt with in Section [tI 



For the time being, we suppose that we have the solutions s^^\nN, ■ ■ ■ ,ni) and 
A'^\n]\[, . . . ,ni) to examine gauge invariants. We also assume that the spectral 
degeneracies are absent, i.e., p2, . . . ,Pn are odd, as explained in the previous section. 
The permutation matrix &^^^ is precisely determined by s^^\nN, ■ ■ ■ ,ni): 

N 

i=i 
Note that the correspondence between {n^, . . . , ni) and rriN^nN, ■ ■ ■ , ni) are one-to-one. 
When we employ the "mAr"-representation, instead of the {un, ■ ■ ■ ,ni) representation, 
the matrix elements of G^^' have simpler forms. Because the change induced by a cycle 
C is represented by the shift of m^v by (i at (mod 2^), we have 

{® }m',niN = ^m',mN+dN mod 2^- (48) 



'JV> 



When d^ is odd, or equivalently when pi is odd (see (41)), mAr(mod 2 ) returns to 
the initial value only after 2^ repetitions of evolution along C. Hence the permutation 
matrix & describes a cycle. On the other hand, if d]\f is even, 171^ (mod 2^) returns 
with smaller number of repetitions and the period depends on dN- 
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6. Eigenspace anholonomy 

We examine tlie eigenspace anholonomy for U^'^\\) along the adiabatic cycle C. We 
find the holonomy matrix M^^\C) of U^'^\\) and the geometric phase 7(C) (29) for 



the adiabatic cycle C based upon our gauge theoretic formalism described in Section [3] 
For the sake of simplicity, we focus on the case that all pj's are odd in this 
section, which guarantees that the spectral degeneracy is absent. Each eigenvalue and 
eigenspace returns to their initial ones only after 2^ cycles along C. In other words, 
this assumption implies the presence of the one-to-one correspondence between quantum 
numbers {riiy, . . . , rii) and m^{njq, . . . , Ui). Hence we will often abbreviate (tzat, . . . , ni) 
as rriM in the following. Also, for the sake of brevity, we will drop the index C. 

As shown in Section |3| M^^^^ can be decomposed into a product of the permutation 
matrix (S^^-* and a diagonal unitary matrix, i.e., 

{mW}„^,^, = {6(^)}^^,^,a(^)(m^), (49) 

where a^'^\mN) is a unimodular complex number. Because &^^^ has been already 
obtained in the previous section, our task here is to determine the phase factor a^^^rriN). 



Let us obtain a recursion relation of M^^^ ((59 ) below) using the gauge connections. 
The eigenvectors of f/^(A) form a non-Abelian gauge connection 

{A^^'^X^m'mr, = K^; M WniN; A)] , (50) 



and the diagonal connection {A^ >{X)} 



m',m,fi 



= Sr, 



,m]v 



{AW(A)} 



m'i^,mi^i 



which are 



2^ X 2^ Hermite matrices. We now obtain a recursion relation for A^'^^X). For A^ = 1 
(4a), (Ab) and the definition of the gauge connection imply 

2 



A«(A) 



P 



4 



Y, and A°(i)(A)=0. 



(51) 



Namely, for A^ = 1, each eigenvector satisfies the parallel transport condition within 
each eigenspace. For A^ > 1, the Leibniz rule suggests a decomposition ^(^^(A) = 



iW 



(N), 



A[^^^(A) + <'^(A),where 



{4"Ha)} 



m',mr^ 



^{AiQ'iX),p^)},d,Q'{X)5, 



^jV-l'™'iV-l 



m'j^,mM 



(52) 
(53) 



n'^(e'(A),p;v)|n^(e(A),p^)){A(^-i)(A)}^^_^,^^_^ 

Here e(A) = O^^-^^mN-i; X) and e'(A) = e^^~^'^{m'^_^; X) are introduced for 
abbreviations. We also have a similar recursion relation for "diagonal" gauge connection 
A^^^\X). Because we have chosen the gauge ^'^'^^^(A) = for A^ = 1 (see, (51)), it is 
straightforward to obtain A^^^\X) = for an arbitrary A^ from the recursion relations. 
We obtain an explicit expression of A^ (A) from the eigenangle of A^-qubit 
system (38) and the gauge connection of the single-qubit system (51): 

(2 



iW 



<'W 



-PN)d 

2N+1 



N-1 



r® 1(^-1). 



iW 



Next we examine A\^ (A) in our model. From (4a), (46) and (53), we obtain 



iW 



A^'W 



1 + Y 



A 



(N-l) 



(A) 



1-Y 



A 



(N-l) 



(A), 



(54) 



(55) 
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and 



J^ 



m',mff 



Accordingly we have 



= 171^6. 



mj^,mM- 



2-'vy®j(f-i' 
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(56) 
(57) 

(58) 



Now it is straightforward to prove that A^ (A) and A^^\\) are independent of A. First, 
we assume the independence of A^^~^\X) from A for A^ > 1. This assumption imphes 
that the A^ (A) is also independent of A, from (58). Because A\^ (A) (|54[) is also 



independent of A, we conclude the independence of ^(^^(A) from A. Second, A^^\X) is 
independent of A (see, (51)). Hence we conclude A\^ (A) and A^^\X) ioi N > 1 are 
independent of A. 

In addition, it is straightforward to see [A]^ '(A), A\^ (X)] = 0. Hence the anti-path 
ordered product is simplified as M^^'^ = q-^'^'^'^h Wq-^'^'^-^-l W, This implies a recursion 
relation for the holonomy matrices: 

^ + ^ - ^>(2-pjv)2-^J^'^-'>j\^(iV-l)g-i7r(2-pjv)2-^(J^^-l>+dAr_i)A 



M(^) 



2 

1-Y 

H — (X> |e 



i^(2-p^)2-^j(,'^-i)^(7V-l)gi^(2-pjv)2-^(j(f-''+rfiv-i)A _ 



(59) 



From this recursion relation, we obtain the recursion relation for the matrix element of 

Mft 



{M<'')}„, 



e 



(._.,)(-l)^("-P-^"-^'){M(^-i)}^, 



jV_ii'"iV-l 



(60) 



which is shown in Appendix A 



We now obtain a^'^\mj^), which is a constituent of M^'^^ (see (49)). For A^ = 1. it 

(61) 



is easy to see a^^\ni) = (— 1)''('^1'Pi). From (60), we find 



m'j^,mM 



{G^^'^U>^A-iY V 



^W (Ar_l) 



rriN-i, 



which implies a recursion relation for a^^^: 



a 



Wi 



(m^) = (-l)^^"'(-^)a(^-i)(m;v-i). 



(62) 



Hence we have 

a(^Hnjv,...,ni) = (-i)!:^^!-*'^".,-,-!)^ (63) 

where we set r^^\ni) = r{ni,pi). Note that an explicit expression of r^''\nk, ■ ■ ■ ,ni) 
depends on p/s. 



We examine Manini-Pistolesi's gauge invariant 7*^^^ From (24), we obtain 
1 1 



^17 



(iV) 



njv=0 ni=0 



(64) 
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In particular, we have already obtained e^''' = —1 for A^ = 1. For A^ > 1, (62) implies 
ni-_n aW(n;v, . . . , rii) = ni,=o(-l)''"'^""'-'"^^ due to {a(^-i)(n^_i, ...,n,)}'=' 



injv -, 

Thus we have 



= 1. 



3i7(-) = (_r 



where 



AN) 



r,.,— n r,,— n 



(65) 



(66) 



Pnt 



n^=0 ni=0 

Recall that v\' = pi for A^ = 1. From (466), we find ^„j^=o '"''^''(^^' • • • ''"'i) ~ 
■^~^^{n^-\^ . . . , n\) for A^ > 1, because [j/2] + [(1 + j)/2] = j holds for an arbitrary 
integer j. It leads us to a recursion relation v^^^ = pnu^'^^^\ whose solution is 

i/(^) = dN. (67) 

Since we consider the case that no spectral degeneracy exists implying d^ is odd, we 
find 

-1 (68) 



(N) 



for A^ > 0. 



7. Examples 

We examine several examples in this section. To complete the characterization of the 
eigenangle and eigenspace anholonomies of [/'^^•'(A), we need to obtain the explicit 
expressions of s^^h 



UN, 



,ni 



and r*^^)(nAr, . . . ,ni), as explained in Section [sj This 
requires to solve (46a) and (466), whose solution precisely depends on the set of integer 



,N 



parameters {pj} -^i- First, we consider the simplest case pj = 1 for all j. Because we 
obtain the subsequent cases through the modifications of the simplest one, the study 
of the first case offers the basis of the following analysis. Second, we replace pi of the 
simplest case with an even integer. It is shown that the permutation matrix &^^\C) 
of the second case consists of two cycles. Third, we replace pj in the simplest case 
with an odd integer. These examples indicates that there are various types of &^^\C) 
associated with the choices of {pj} •^^• 

7.1. Pj = I for all j : the simplest case 

We here examine the simplest case of pj = 1 for all j (see figures [2] and [s]) , which can 



be reduced to a quantum map under rank-1 perturbation [TJ IB]- From (40), we have 



(69) 



N 

niNiuN, . . . , ni) = ^ 2^~^nj, 
i=i 
which means that n^v, . . . ,ni form coefficients of the binary expansion of the principal 
quantum number ttiat. This provides a direct relation between m^ and {tim., • • • , ""-i). It 
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,i^y){y\ 



A{-l+Pim{y\ 



Figure 2. The quantum circuit U'^"^\\) for three qubits with P2 — P3 - 
orthogonal to \y) and normalized. We examine the cases pi = 1 and 2 in 



7.2 respectively. The first examples shown in 



1. \y) is 
[TTJand 
7.3| correspond to the case pi = 3. 



The convention for the controlled qubit is explained in the caption of figure IT] 




Figure 3. The eigenangles 6'(^)(n3, ^2, rii; A) (38) of U^^^{X) with p2 



P3 



1 (see 



figure pi. The labels of the vertical axis indicate the quantum numbers (?i3, ?i2: "-i) 
of the eigenstates at A = 0. Three families of lines correspond to pi = 1 (bold), 
Pi — 2 (dashed) and pi = 3 (dash-dotted). Note that pi determines the slope of 
(^ {^3, n2,ni; X) . Hence different pi induces different itineraries of quantum numbers 
along adiabatic cycles. In § |7.1[ we examine the simplest case pi = 1. When pi is an 



even integer, 6(C) is decomposed into multiple cycles (see (74) and ( [75| ) in § 7.2). On 
the other hand, odd pi ensures that &{C) describes a cycle (see ([83]) i 
Pi = 3 is assumed). 



7.3 where 



is straightforward to solve (46a) and (466) 



JN) 
'AT 



sm 



(nAr,...,ni) 




for ni\i_i ■ 
otherwise 



r)y ^(riAT, ...,ni) = n^ ■ ■ ■ ni, 



where we define the product of quantum numbers as 

for k > j 
1 otherwise 



rife ■ ■ ■ Uj 




ni 



(70 a) 
{70b) 



(71) 
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throughout this paper. 

The principal quantum number rrii^ increases by one after parametric evolution 
along C once due to d^ = 1: tjin- = H- 1 ^-)■ . . . 2^~^ — 1 h-)- ^-)■ . . ., whose period is 
2^. On the other hand, the itinerary of tin, . . . ,ni oi N = 3, for example, is 

000 ^ 001 ^ 010 ^ Oil ^ 100 ^ 101 ^ 110 h^ 111 ^ 000 h^ . . . . 

7.2. The simplest case perturbed with pi ^ 1: pi is even 

We examine the case of pj = 1 for j > 1 and pi is an even integer, which corresponds 
to the simplest case perturbed with a certain even pi ^ 1 (see figures [2] and [s]) . Note 
that the mapping between n^v, . . . , ni and m^v is independent of pi, while the spectrum 
depends on pi. 

9^^\nN,...,ni) = ^ (m7v(njv,...,ni)+pi— j . (72) 

Hence as A increases by every 2tt, the principal quantum number rriN also increases 
by pi. It means that rriM mod 2^ returns to the initial point after 2^~^ cycles of C. 
Furthermore, a suitable choice of pi make the period even shorter. This is the crucial 
difference from the case that pi is odd, as discussed later. 

For the sake of simplicity, we assume pi = 2 so as to find the explicit expressions 
of s^^^ and r^^\ For A^ = 1, we have Si (^i) = ni, r^^\ni) = 1. For N >2, one finds 

(TV). s I 'W ioT riN^i ■ ■ ■ n2 = 1 

s}^'{nN,...,ni) = <; _ , (73a) 

otherwise 




r^^\nN,...,ni) = nN---n2. {73b) 

Hence the first qubit is decoupled from others. The rest of the qubits are equivalent to 
the simplest case with A^ — 1 qubits. 

For pi = 2, the permutation matrix &^^\C) forms two sorts of cycles consisting of 
even or odd mjy's. For A^ = 3, the case of even rajy reads 

000 ^^ 010^ 100^110^000^ ..., (74) 
while for odd mjv's 

001 h^Oll ^ 101 h^ 111 h^OOl f^ .... (75) 

We remark that ni remains unchanged along these itineraries. Namely, the initial qubit 
is repeatedly recovered during the evolution. 

7.3. The simplest case perturbed with odd pj 

Here we consider the case that pj = p6jj + (1 — Sjj), where J is an integer (0 < J < A^) 
and p = 1 + 2^^ {K is a positive integer). In other words, we obtain this model by 
replacing pj of the simplest case, where we set pj = 1 for all j, with p. We examine 
how such a tiny change affect the eigenangle and eigenspace anholonomies. When this 
new model consists of smaller number qubits, i.e., A^ < J, it becomes equivalent to the 
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simplest model. The effect of the replacement of pj appears only when the system size 
is large, i.e., N > J. 
First, we note that 

N J-1 

m^inN, ...,ni) = Y, '^''\ + Y. '^^^'~^^v (76) 

where the second term in the rhs is defined as if J = 1. We also note that the spectral 
degeneracy is absent since p is odd 

We show explicit expressions of s)q {um-, ■ ■ ■ , rii) and r^^\nN, . . . , ni), which govern 
the anholonomy in the quantum number. Because this mode l agr ees with the simplest 
case for N < J, s]^ (^Af, . . . ,ni) and r}^ {n^, ■ ■ ■ ,ni) satisfy (70a) and {70b) for N < J. 
The effect of the pj ^ 1, which we may call an "impurity" of the simplest model, sets 
in at A^ = J. For J < N < J + K, one finds 

(TV) . \ i^ ^°^ ^^-1 ■■■^1 = 1 .^7 ^ 

s}^'{nN,...,ni) = < , , (77 a) 

I riN otherwise 

r'~''\nN, ...,n,) = {nN---nj + 2^-(^-^)-i)n,,_i ---m. {77b) 

Hence the effect of the impurity pj on r^^^ is largest at A^ = J, and becomes smaller as 
A^ increases for J < N < J + K. In other words, r^^^^s of K qubits are directly affected 
by the impurity. The J + K-th qubit is in another regime: 

(j+K) , ^ f ^^J+^ for ^J+K-1 ■ ■ ■ nj = and nj_x ...n^ = \ 

s^j^j^\nj+K,---,nx) = \ . . (78) 

I "^j^-K otherwise 

The influence of the qubit with index j (j = J, . . . , J + K) on the qubits with larger 
index can be described through a quantum number 

tiJ-K) {nj+K, ...,nj)= nj+K + nJ^nj+K-i ■■■nj, (79) 

which is either or 1. It seems that the impurity introduces an effective {K + 1)- 
body interaction among the quantum numbers {tij+k, ■ ■ ■ ^nj). This determines r^^-* 
for iV > J + ir as 

r^^\nN, . . . , ni) = nTv . . . nj+K+it^'^'^^nj^i ■ ■ ■ m, (80) 

where we omit the argument of t^'^'^\nj^K, • • • ? nj). On the other hand, for N > J + K, 
we have 



, um for riTv-i ■ ■ -n j+K+it^'^'^^'n j.i • ■ ■ rii = 1 
'^'(n^,---,m) = <l "" ^^ ""/ '^''^" . (81) 

TiTv otherwise 



SN)i 



We consider two typical cases. The first example is the case of J = 1 and K = 1 (see 
figures [2] and |3]) , where the principal quantum number rriN is unaffected by the impurity 
due to J = 1 (see, (76)). As shown above, a cycle C corresponds to the increment of 
iriM hy (In = Pi = "i- For example, with A^ = 3, we have 

0h^3^6^1^4^7^2^5^0^.... (82) 
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f/(i)(A) 


— = — 


Z 




e'^|v>(>i 



[/(2)(A) 



T 



[/(i)(A) -{[/(i)(A)}2 




f/(i)(A) -{t/(i)(A)} 



Figure 4. The quantum circuits ^/'^•'(A) (iV = 1,2, and 3) witli pi — ps — 1 and 
P2 = 3. Tlie cases of A^ = 2 and 3 are ratlier complicated due to tlie recursive 
structure with p2 ^ 1. This is the second example shown in § |7.3| The convention for 
the controlled qubit is explained in the caption of figure [l| 



In terms of njS, the itinerary is 

000 ^ Oil ^ 110 ^ 001 ^ 100 ^ 111 ^ 010 h^ 101 ^ 000 h^ . . . . (83) 

This itinerary can be understood through the simple correspondence between n^, . . . , ni 
and mi\i, as shown in (76). Namely, njsf, . . . ,ni are coefficients of the binary expansion 
of rriN- 

The second example is the case of J = 2 and K = 1 (see figure |4] and Isl) , where 
the effect of pj 7^ 1 appears not only in the eigenangle but also in niN^nN, • • • , ni) (see, 



(76)). On the other hand, the itinerary of rriN induced by the repetition of the cycle C 
remains intact to be equivalent to the case J = 1 and K = 1 discussed above; i.e., a 
cycle C corresponds to the increment of rriN by 3. We depict the itinerary of quantum 
numbers Hj induced by the cycle C. For A^ = 3, we have m^ = An^ + 2n2 + 3ni and 

000 ^-> 001 y-> 110 H^ 111 H-> 100 ^ 101 ^ 010 h^ Oil ^ 000 ^ . . . . (84) 



8. Complexity of the anholonomy 

The repetition of the adiabatic cycle C in U^^\\) generates an itinerary of the quantum 
numbers {un, • • • , ni). We point out that the itinerary involves an NP-complete problem 
in this section. 

The examples in the previous sections tells us that the expressions of s]^ s and 
r'^^^'s are rather simple when pj = 1 for all j. In contrast, when we make any pj 
different from 1, the expressions of Sj^ s and r'^^^'s become slightly complicated. One 
may expect that the more p/s differ from 1, the more complicated they become, so 
does the itinerary of {n^, . . . ,ni) induced by the adiabatic cycle C. However, this is 
not true because the itinerary of ttiat remains simple irrespective of variation of p/s; 
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Figure 5. The eigenangles 0''^^{n3,n2,ni; X) (38 1 of f7^^^(A) depicted in figureUl The 
labels of the vertical axis indicate the quantum numbers (713, n2,ni) of the eigenstates 
at A = 0. The order of the quantum numbers at the initial point A = is different 
from that of figure [3J Namely the initial quantum numbers are scrambled by p2 7^ 1- 
The corresponding itinerary of quantum numbers is shown in (84 1. 



TTiN increases only by d^ for a cycle C. How do we justify this? We will show that 
the cost to obtain n^ ... ni for a given rrijy is generally equivalent to that of finding the 
solution of the subset-sum problem, which is a NP-complete problem. We also explain 
the relationship between their equivalence and the adiabatic quantum computation |T3] 
along an eigenangle [H]. For the sake of simplicity, we focus on the case that all pj is 
odd. Then, the spectral degeneracy is absent, and the shortest period of the itinerary 
is 2^ in the unit of the cycle C. 



First, we explain our task to find {njy . . . rii) for a given rriN. From (40), this is to 
solve the following equation 



N 



rriN 



^Wj-rij, 



^5) 



Wj = 



where 

N \ 

n P^ 2^"i (86) 

\k=j+l J 

for 1 < j < A^, and wn = 2^~\ 

Second, we explain the subset-sum problem [22j. A subset-sum problem has two 
parameters S and t. S' is a finite set of positive integers, t is a positive integer, and is 
called a target. The problem is to show whether there exists a set S" C S* such that 
t = Xlses" ^- -^^^ °^^ purpose, we need to obtain S'. 

Now we explain the equivalence between our problem and the subset-sum. In our 
case, let S be {wj}.^, which is a set of positive integers. Finding {n^- . . .ni) that satisfy 
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(85) is essentially equivalent to the problem to find a set of positive integers S" C S* 
such that 171^ = ^sg5' s. This justifies the equivalence. 

This explanation has a subtle point. During a course of the repetition of the cycle 
C, m^ can be determined only up to modulo 2^. However, this does not break the 
equivalence in the difficulty of the solving the problems 

Based upon this equivalence, we examine whether the evolution along the adiabatic 
increment of A in U'^^\X) offers an efficient way to solve the subset-sum problem. It 
turns out that the adiabatic approach is inefficient. First, when A is varied slow enough, 
we can prepare irip^ to be an arbitrary positive integer. Next, appropriate measurements 
of qubits ascertain n^v, . . . , ni. To make the nonadiabatic error small enough, there is a 
lower bound of the running time for quantum evolution. This is essentially determined 



by the gap A between the eigenangles. From the exact expression of the eigenangle (38), 
we obtain A oc 2~^ for large N . Hence the lower bound of the running time must be 
exponentially large, i.e., the adiabatic approach is inefficient. 

We remark on the inefficiency to obtain tj'^^\X) from quantum gates. Our 



construction described in ([33J), of tj^^\X) is generally inefficient, because (33) requires 
exponentially many quantum gates for general {pj}°°^q- On the other hand, if we impose 
Pj = 1 for most j except for infinitely many j, the number of quantum gates required 
to construct U^^\\) can be a polynomial of A^. However, we are not certain whether 
such a U^^\\) involves any NP-complete problem. 

9. Discussion 

We here discuss the relationship between our result and the recent work on the 
topological characterization of periodically driven systems [T3]. The work in [IS] offers 
an integer 

u{C) ^^i Tr[{f/(A)}-ia,f/(A)]dA, (87) 

2vri Jc 

for a family of Floquet operator f/(A) in a closed path C in the parameter space of A. 
Because z^(C) can be regarded as a gauge invariant for the adiabatic cycle C, it is worth 
comparing it with our gauge invariants ©(C) and 7(C). 

First, we show that uiC) is indeed an integer. Let U{\) denote a Floquet operator 
of a periodically driven system. We assume that f/(A) is periodic in A and its shortest 
period is 27r, which guarantees its spectrum je'^"^'^-'} _ is also periodic in A, where 
A^ is the dimension of the Hilbert space and On{\) is an eigenangle. Such a spectral 
periodicity implies 

^sH(0)=e„(27r) + 27rr(n), (88) 

where s{n) describes a permutation over quantum numbers {0, . . . , A^ — 1}, and r{n) is 
an integer. The latter integer has an integral expression 

2.r{n) = r ^dA. (89) 
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The derivative of eigenangle is d9n{X)/dX = —i{n{X)\{U{X)}~^{dxU{X)}\n{X)), where 
|n(A)) is a normahzed eigenvector of U{X) corresponding to the eigenangle 6'„(A) [23] . 
Hence we obtain 

/■27r 

2Txr{n) = -i {n{X)\{U{X)}-^{dxU{X)}\n{\))dX. (90) 

A sum rule on 27rr(n) has a representation- independent expression: 

27rVr(n) = T/ Tr {U{X)}-\dxU{X)} dA. (91) 

n=o 1 io L J 

We remark that this sum rule is the key to prove the presence of the eigenvalue and 



eigenspace anholonomies in the quantum map under a rank-1 perturbation [HI [7]. (87) 



and (91) imply 

Af-l 

^(C)^^r(n), (92) 

n=0 

which is an integer. This argument suggests that ^{C) characterizes the winding of 
eigenangles, or equivalently, quasienergies. 

We now compare ^{C) with the permutation matrix &{C). We examine the 
classification of a closed path C for a family of the quantum circuit on a qubit, using 
^{C). We recall that ©(C) is shown to classify C into two classes in Sec. [2J This is 



obtained by an inspection of Cp, which corresponds to the family of unitary u(X,p) (pi) 



for < A < 2n. In particular, even and odd p correspond to (12) and (14), respectively. 
As for z/(C), it is straightforward to see i^{Cp) = p. The stability of the topological 
quantity z^(C) against a small deformation of closed path from Cp implies that i^{C) of 
a closed path C takes an arbitrary integer. Hence z^(C) classifies the closed paths into 
an infinitely many classes. Thus z^(C) exhibits a detailed structure than &{C) as for 
the closed paths for the space of single-qubit quantum circuits. 

This conclusion does not hold for the quantum circuits with multiple qubits. To see 
this, we examine the quantum circuits U'--^\X) introduced in Sec. |4| We have already 
obtained the corresponding integer u^^\C) in (67) to evaluate the geometric phase 



■y(^)(^C). The examples shown in Sec. 7.3 tells us that &^^\C) can take various values 
even for a certain given u^^\C). Thus we conclude that the role of ©'-^-'(C) and u^^^C) 
are independent in classifying the families of quantum circuits. 

In order to provide another comparison of &{C) and z^(C), we show an example in 
which C contains the crossing of eigenvalues. Note that we have excluded such cases so 
far in this paper. We will show that &{C) is sensitive to the eigenvalue crossing, while 
z/(C) is not pK]. Let us examine the following quantum circuit 

UYiX,p) = e'(P"i)^(^-l^><^l)+'^l^><^lf, (93) 



which is obtained by replacing Z with Y = 1 — 2\y){y\ in u{X,p) (|2j). These quantum 
circuits are periodic in A. Let Cy denote the closed path in quantum circuits specified 
by uy{X,p) with < A < 27r. The crucial point is that Cy involves an eigenvalue 
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crossing, because the eigenvalues of ^^(A,^) degenerate at A = 0. Also, the eigenvector 



of uy{X,p) is independent of A. Hence &{Cy) is the identity matrix (12). We explain 
that &{C) in the vicinity of C = Cy is generically different from the identity matrix. 
Let us introduce another unitary matrix Y', which satisfies [Y', Y] ^ 0. When we 
replace Y in uy{^,p) with Y', the corresponding family of quantum circuit exhibits the 
eigenspace and eigenangle anholonomies, according to the analysis in [6]. Hence &{C) 



is the permutation matrix whose cycle is 2 (14). This conclusion holds even when the 
difference between Y and Y' is arbitrary small. Thus, &{C) is sensitive to C in the 
vicinity of Cy. In contrast, z^(C) = p is independent of Y in uy{X,p), as long as Y is 
unitary. In this sense, i^(C) is stable against the choice of C, even when C involves the 
crossing of eigenvalues. 

We close this discussion with a comment on the gauge invariants obtained here. All 
gauge invariants are determined by the integers s{n; C) and r{n\ C), which are defined 



in (43) (see also, (88)), as for the adiabatic cycle C of quantum circuits U^^\X) (33). 



The permutation matrix ©^^^ contains all s(n)'s. On the other hand, i^{C) (92) is the 



whole sum of r(r2)'s. We find that i^{C) determines the geometric phase 7(C) through 



(65). It remains to be clarified whether the intimate relationship between ^{C) and 
7(C) holds in general. Also, it is worth to clarify whether other combinations of r(n)'s 
give us any useful insights on the anholonomies. 

10. Summary 

In this paper, we have identified the gauge invariants &{C) and 7(C) that are associated 
with the eigenvalue and eigenspace anholonomies for the closed path C of a family of 
unitary operators. The unified theory of quantum anholonomy has been revisited to 
clarify that 7(C) is the Berry phase for the A^-repetition of closed path C^ , where A^ 
is the dimension of the relevant Hilbert space. By using a family of quantum circuits 
that are recursively constructed, these gauge invariants have been analyzed in detail. 
It has been shown that a generic family of the quantum circuits is associated with an 
NP-complete problem. The relationship between our gauge invariants and Kitagawa et 
a/.'s topological integer i'{C) has been also discussed. 
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Appendix A. A proof of (60) 



recursion relation (59) of M^^^ 



We will obtain the recursion relation (60) of the matrix element of M^^^ from the 



(AT) 



AN) 



First, to simplify (59), we need to evaluate e'^'^n M'-^'^e ^"""^d . Because Jf, (57) 



is diagonal and M^^'> is the product of a permutation and a diagonal unitary matrices 
(see, ([49])), we have 

{[J^''\ M(^)]| = {M(^)}^.^,^, [m^(s(^)(n^, . . . ,m)) - m^{n^, ..., n,)] .(A.l) 



m'j^,mM 



Furthermore, using the recursion relation (|45|) for rrijv, we obtain 



{M(^^U--iv(^^-2V^)), (A2) 

where r^^^riN, . . . ,ni) is abbreviated as A'^\ Now it is straightforward to show 



m',mM 



giaj(-)^(7V)g-iaj(-' 



m'^,mN 






Hence, we obtain from (59) 



m'j^,mM 



i(V2)(2-p^)r('V-i)y| (m(^-^)) 



n^,nf^ 



"lW-1'"^iV-l 



When r*^^ -^^ is even, the first factor in rhs of (A. 4) is written as 



-i7r(2-pjv)r(^-i)y/2\ 



J n' " 



,raiV 



On the other hand, when r*-^ ^^ is odd, we have 



-i7r(2-pjv)r(^-l)y/2l 



^r! 






\nN ( _T\[{2-PN)r^^-''HnN-i,-,ni)-l]/2 



We then obtain, using the definition of r{n,p) ^ 

-i7r(2-pjv)r(^-i)y/2 



n^,nff 






Hence (60) is proved. 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 
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